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ABSTRACT 

We report on the present stage of SN 1987 A as observed by the Chandra X-Ray Observatory. We reanalyze 
published Chandra observations and add three more epochs of Chandra data to get a consistent picture of the 
evolution of the X-ray fluxes in several energy bands. We discuss the implications of several calibration issues for 
Chandra data. Using the most recent Chandra calibration files, we find that the 0. 5-2.0 keV band fluxes of SN 
1987A have increased by ~6 x 1 0 1 ' erg s _1 cm -2 per year since 2009. This is in contrast with our previous result 
that the 0. 5-2.0 keV light curve showed a sudden flattening in 2009. Based on our new analysis, we conclude that 
the forward shock is still in full interaction with the equatorial ring. 

Key words: ISM: supernova remnants - radiation mechanisms: thermal - supernovae: individual (SN 1987 A) - 
X-rays: individual (SN 1987 A) 
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1. INTRODUCTION 

Stellar winds of massive stars leave their imprints on the 
surrounding medium before they explode. After the explosion, 
the outer shocks of the supernova probe the circumstellar 
medium, thereby mapping the late stages of the stars evolution. 

Because of its youth, the remnant of SN 1987A provides a 
unique opportunity to investigate these late stages. The case of 
SN 1987 A is especially interesting, as it exhibits the enigmatic 
three ring structure visible in optical images (Crotts et al. 1989). 
This structure probably originates from the interaction of a wind 
from the slow red supergiant phase with the faster wind from the 
blue supergiant phase (Luo & McCray 1991; Wang & Mazzali 
1 992). The morphology of this structure might provide clues to a 
binary-merger scenario for the progenitor (Blondin & Lundqvist 
1993; Morris & Podsiadlowski 2007, 2009). However, this three 
ring structure might also result from mass loss from a fast- 
rotating star (Chita et al. 2008). Currently, the explosion is 
sweeping up the inner equatorial ring that was formed by the 
late stages of the star’s evolution. 

As the inner equatorial ring is about a hundred times denser 
than its surroundings (e.g., Mattila et al. 2010), the interac- 
tion of the explosion with this structure has a direct impact on 
the evolution and observational properties of the remnant. In 
1995, the outer ejecta began to light up this structure. By 2000 
October, 12 spots had lit up at visible wavelengths (Sugerman 
et al. 2002). Infrared observations indicate that the full interac- 
tion with the ring started around mid-2002 (Dwek et al. 2010). 
A year and a half later, the expansion in X-rays slowed down 
dramatically (Racusin et al. 2009) while, at the same time, the 
soft X-ray light curve turned up (Park et al. 2005). 

Based on past Chandra observations, Zhekov et al. (2005) 
established that the X-rays are emitted from a flat spatial struc- 
ture. A model consisting of three components has been proposed 


for interpreting the high-resolution X-ray spectroscopy from 
SN 1987A (Zhekov et al. 2010; Dewey et al. 2012). The first 
component is emitted by material in the dense ring that has been 
heated by a slow shock (500 km s _1 ). This component has a low 
electron temperature (~0.3-0.5 keV). The second component 
is emitted by the shock-heated circumstellar medium located 
above and below the equatorial ring. These shocks have not 
been slowed down as much and therefore the plasma is hotter. 
Also, the line emission of this component has to be broad (thou- 
sands of km s _1 ). This allowed Dewey et al. (2012) to estimate 
the current contribution of this component to be of the order 
of ~20% of the total emission in the 0. 5-2.0 keV band, based 
on high-resolution High Energy Transmission Grating (HETG) 
spectra and Reflection Grating Spectrometer spectra. Further in- 
vestigation of the X-ray spectra of SN 1987A shows that there is 
need for a third component with a higher (~3 keV) temperature, 
but without broad-line emission. Zhekov et al. (2006) proposed 
that this is shocked plasma, reheated by a reflected shock caused 
by the interaction with the equatorial ring. Alternatively, Dewey 
et al. (2012) proposed that this high temperature component 
arises from a fast shock, traveling through the uniform equato- 
rial ring material (in contrast to the low temperature component, 
which they consider being emitted by shocked denser clumps in 
the equatorial ring). 

The Chandra X-Ray Observatory has observed SN 1987A 
at a regular pace of about two observations per year, resulting 
in a denser coverage of the evolution of the remnant than with 
any other X-ray observatory. Park et al. (2011) summarized 
the observations through 2010. In the 0. 5-2.0 keV light curve, 
they found a break (i.e., the pace with which the flux increases 
suddenly slowed down) around 8200 days after the explosion 
(~Y2009). Based on this result, they concluded that the shocks 
responsible for the soft X-ray emission must have recently 
encountered less dense material. 
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ACIS-S3 Observations 
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Instrument 

Frame Time 
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3.0-8.0 keV Flux a 
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(s) 

( 10 13 erg s -1 cm' 2 ) 

(10~ 13 erg s _1 cm -2 ) 
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3.2 

2.79!° c j. 1 1 1 3 (L18) 
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0.4 

17 - 71 -d 6 M l 1 ' 08 ) 
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0.4 

21 . 827 °o 7 8 6 5 ( 1 . 11 ) 

2-537*0.25 (°- 94 ) 
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6914 

ACIS-S3 

0.4 

26.727 ° 0 7 9 6 9 (1.13) 

3.237°o 1 2 8 s (0.96) 
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7094 
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0.4 

3\.6C^\ (1.14) 

3.1 (0.85) 

7636 
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0.4 

37.687°o 4 4°3 (1.18) 

3. 737*o.\*7 G- 89 ) 
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7445 
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42. 277q 29 (1-19) 

3.397°o 1 1 4 4 (0.77) 

9142 and 9806 
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ACIS-S3 

0.2 

46. 377o4 5 6 (1-09) 

3.427°o 2 3 1 6 (0.63) 
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7799 

HETG 

1.1 

50.367 \' 5 7 6 0 C- 06 ) 
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1.1 
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4.867° 4 1 o 3 6 (0.81) 
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8233 

HETG 

1.1 

60.787 2 3 4 6 5 6 (1.07) 

5.837o. 7 2 (0.89) 

12125 and 12126 and 11090 

8434 

HETG 

1.0 

63.507 2 3 6 4 2 8 (1.08) 

5.957*o.4i (0-90) 

13131 and 11091 

8618 

HETG 

1.0 

66 1 5 -1.52 (!■ 0 8) 

6.087° 0 23 o (0.85) 

12539 

8796 

HETG 
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Notes. Observation log of the Chandra S3 observations of SN 1987 A. 

a Numbers between parentheses are the ratios of the pileup-corrected fluxes and the non-pileup-corrected fluxes. 


However, the 2012 May revision to the ACIS calibration 
(CALDB 4.4.1 0) : introduced a significant change to the time- 
dependent model for the transmission of the optical blocking 
filter. The new model predicts lower transmission than the 
previous model (a 15% change in the optical depth in 2012), 
reflecting an apparent increase in the rate of deposition on the 
filter. 

In this paper, we reconsider the light curve break result (Park 
et al. 2011) by reanalyzing previously published epochs using 
the newest calibration and by adding three recent observations 
(2011 March and September, 2012 April). We show that the 
break seen by Park et al. (2011) is due primarily to the 
(unmodeled) increased rate of deposition on the optical blocking 
filter and not to a different phase in the evolution of the remnant. 

2. DATA AND RESULTS 

Table 1 lists the 23 monitoring observations used for this 

study. We chose not to use the observations earlier in the mis- 

sion (Obslds 1387 and 122 at 4608 and 4711 days after the 

explosion), as these observations were taken at a focal plane 

temperature of —110° C and no charge transfer inefficiency 

adjustments are available for observations taken at this temper- 

ature. 9 However, we report the spatial distribution of the X-ray 

data at these epochs in Section 2. 1 . Since the remnant has thus 


8 http://cxc.harvard.edu/caldb/downloads/Release_notes/CALDB. 
v4.4.10.html 

9 http://cxc.harvard.edU/ciao4.4/why/cti.html 


far brightened by a factor of ~34 during the Chandra mis- 
sion, the observing configuration was changed several times to 
mitigate photon pileup effects. 10 First, this was done by read- 
ing out an increasingly smaller subarray of the ACIS-S3 chip. 
All observations from day 7799 onward used the AC1S-HETG 
configuration, 11 which produces both a zeroth-order image and 
a dispersed spectrum. The dispersed spectrum provides an in- 
valuable cross-check for our ACIS results. From day 8434, the 
observation was moved to the bottom of the chip, to decrease 
the frame time to 1 s. The time of this transition coincides with 
the accelerated pace of the increase of the contamination on the 
filter. For the data reduction, we largely follow the methods that 
were used before on Chandra data of SN 1987A (Burrows et al. 
2000; Racusin et al. 2009; Park et al. 2011). We filtered for 
background flares, and ignored the flags set by the Chandra af- 
terglow filter (status=000000000000xxxx0000000000000000) 
as this filter tends to filter real events for bright sources (e.g., 
Townsley et al. 2003). 

2.1. Imaging 

With an angular diameter of 1" 5, SN 1987A is barely resolved 
in raw Chandra ACIS images. To extract as much spatial 
information from the observations as possible, we adopt the 
same approach as in previous studies (Burrows et al. 2000; 
Racusin et al. 2009; Park et al. 2011). We remove the pixel 


10 http://cxc.harvard.edu/ciao/ahelp/acis_pileup.html 

11 Use of the HETG is an excellent pileup mitigation strategy, since the count 
rate in the zeroth order is reduced by a factor of nine. 
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Figure 1. Estimated intrinsic morphology (corrected for the Chandra PSF) of SN 1987 A at all epochs utilized in this study, proceeding from upper left to lower right. 
The image scaling is linear and the brightness scales with the 0.5-8. 0 flux in the observation. The scaling limits of the two top rows are offset by a factor of three for 
visibility. In the images, north is up, east is left. 


randomization that is usually added by the Chandra software. 
Additionally, we improved the spatial resolution using split- 
pixel events (Mori et al. 2001). We then deconvolve the 
image for the point-spread function of the telescope 12 using 
the Lucy-Richardson iterative deconvolution algorithm (Lucy 
1974; Richardson 1972). 

For the 2011 September and 2012 March observations, we 
analyze only the longest exposure to avoid distortions caused 
by residual alignment errors between exposures. Figure 1 shows 
the deconvolved images of all epochs of Chandra observations 
to date. To characterize these images, we fit the image with a 
model that consists of four lobes and a ring (Racusin et al. 2009). 
Figure 2 shows the best-fit radii of all the images considered in 
this study. 


2.2. Spectra 

Spectra were extracted from each observation within a cir- 
cular aperture (radius 4'.'38) sized to encompass nearly all the 
light from the target in all epochs, using CIAO version 4.4 
(Fruscione et al. 2006) and CALDB version 4.4. 10. Background 
spectra were extracted from each observation within an annular 
region (radius 6"2 to 12'.'4). 

Even though we adjusted the observing modes to minimize 
pileup, PIMMS suggests that a modest pileup is present in many 
of the observations. We first tried to take this into account 
during the fitting process by using the pileup model described 
by Davis (2001). This approach adds to the overall model 
several pileup parameters which then have to be either assigned 
appropriate values or included in the fitting; for details see 
the “Chandra ABC Guide to Pileup.” 13 As an example, the 


12 Using MARX version 4.4.0 (http://space.mit.edu/cxc/marx). 

13 http://cxc.harvard.edu/ciao/download/doc/pileup_abc.pdf. Section 4. 1 : 
Correcting Imaging Observations. 
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Figure 2. Radii of a toroidal model (Section 2. 1) fit to each image in Figure 1 . 


parameter nregions encodes the spatial size of the source in 
units of 3 x 3 pixel cells. Based on the ACIS image, we can 
estimate that nregions should be in the range of 2-5, but we 
cannot determine the exact value a priori. Including nregions 
in the fitting, however, adds a degeneracy between nregions 
and the other model parameters such that the fit and its flux are 
not well constrained. 

An alternative, parameter-free strategy for addressing pileup 
in the analysis of ACIS data has been described by Broos et al. 
(201 1, Appendix A). First, a high-fidelity simulator of the ACIS 
detector (Townsley et al. 2002a), including charge transfer inef- 
ficiency (Townsley et al. 2002b), is used to infer an input pho- 
ton spectrum that would cause ACIS to produce the observed 
spectrum and event rate under the specific configuration of the 
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Figure 3. Original (black) and pileup-corrected spectrum of Obsld 4614 with best-fit model overplotted. Note that the pileup spectrum is substantially harder than the 
pileup-corrected spectrum. 

(A color version of this figure is available in the online journal.) 


observation. That photon spectrum is very flexible; no spec- 
tral model is assumed. The spatial distribution of the simulated 
photons incident on ACIS was produced by the MARX pack- 
age, 14 using the corresponding SN 1987A morphology model 
described in Section 2. 1 . MARX and the CCD simulation were 
configured to match the parameters of the observation, such 
as the target location on the focal plane, frame time, and dead 
area/time. 15 

Second, with the input photon spectrum fixed, the simulation 
is repeated with only a single photon arriving in each frame, 
eliminating pileup. The resulting simulated spectrum, which 
we refer to as “pileup corrected,” is then fit using standard 
methods and calibration files. Figure 3 shows the original and 
pileup-corrected spectrum of the observation that suffered worst 
from pileup (Obsld 4614). We applied this method to all ACIS 
observations reported here. The flux correction attributed to 
pileup in each observation is reported in Table 1. This table 
shows that the correction for the soft band is more smooth than 
the correction of the harder band. This probably indicates that 
our pileup correction method is more stable for the soft band 
than the hard band. 

2.2.1. Light Curve 

We fitted the resulting pileup-corrected spectra with an ab- 
sorbed two-component spectral model, using the X-ray spec- 
tral fitting package XSPEC, version 12.6.0 (Arnaud 1996). The 
spectral model included a component in collisional equilibrium 
ionization (vequil) and a non-equilibrium ionization compo- 
nent (vpshock) using neivers 2.0, updated with atomic data 
for inner shell processes (cf. Badenes et al. 2006). We fixed the 


14 http://space.mit.edu/cxc/marx/ 

15 See discussion of FRACEXPO at 
http://cxc.harvard.edU/ciao4.4/why/acisdeadarea.html. 
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Days after explosion 

Figure 4. Chandra ACIS light curves for the soft (0. 5-2.0 keV, higher curve) 
and hard X-ray (3.0-8.0 keV, lower curve) bands, with pileup correction 
(Section 2.2). XMM-Newton EPIC pn fluxes are shown (brown) for comparison. 
Uncertainties shown are I cr. The XMM-Newton EPIC pn fluxes are displayed 
with 3 a errors. 

(A color version of this figure is available in the online journal.) 

abundances to those found in the high-resolution X-ray Low 
Energy Transmission Grating (LETG) spectra investigated by 
Zhekov et al. (2006) and we fixed the absorption to 2.35 x 
10 21 cnT 2 (Park et al. 2006). Figure 3 shows an example of one 
of the fits. We measured the flux in the 0. 5-2.0 and 3. 0-8.0 keV 
energy bands from the best-fit model for all spectra and esti- 
mated flux uncertainties for a 68% confidence level. Figure 4 
shows the evolution of the fluxes in the low and high energy 
bands for the observations corrected for pileup effects. For this 
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Figure 5. Chandra HETG soft band (0.79-2.1 keV) light curves in the dispersed 
spectra. Note that around day 8200 there is a suggestion of a break. The lines are 
fitted to the data points for similar instrument setups individually. Uncertainties 
shown are 1 a. 
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study, we do not need the high signal to noise that would be re- 
quired for a detailed spectral analysis. Therefore, we decided to 
only correct the longest observation per epoch for pileup even 
if there were more observations available in the same epoch. 
This explains the large error bars of the flux at day 7997, as this 
epoch was divided into five observations. 

For comparison, we overplotted the fluxes of the last six 
epochs of XMM-Newton EPIC pn fluxes (Maggi et al. 2012). 
Note that this instrument has been shown to be extremely stable 
(e.g., Sartore et al. 2012; Plucinsky et al. 2012). 


2.2.2. Chandra Grating Observations 

After the previously described procedures, the 0. 5-2.0 keV 
light curve increases more or less linearly from day 6400 on, but 
there still appears to be a break around day 8200, which does not 
appear in the XMM-Newton EPIC pn light curve (Maggi et al. 
2012). 

Since the pileup corrections (Table 1) after the break are rel- 
atively low (7%-ll%) compared to some earlier observations, 
insufficient pileup correction is an unlikely explanation for the 
break. 

We can confirm this directly using the data taken with the 
ACIS-HETG configuration, since the dispersed spectra in each 
of these observations have no pileup, they provide invaluable 
cross-checks for our pileup-corrected ACIS fluxes. The dis- 
persed spectra were extracted using the standard CIAO tools and 
extraction sizes. The Medium Energy Grating plus and minus 
data are combined and the fluxes in the 0.79-2.1 keV range are 
determined directly from the flux-corrected data, as was done in 
Dewey et al. (2012). In addition to the monitoring observations 
listed in Table 1, we analyzed the deep observations taken in 
2007 (day 7335) and 201 1 (day 8774). These observations have 
been described in Dewey et al. (2012). The light curve obtained 
from the dispersed spectra (Figure 5) shows a break around day 
8200 similar to that seen in Figure 4. Since this break occurs at 
the same time that we moved the aim point on the CCD, we 
cannot disentangle any effects caused by that change from 
those caused by the increased contamination deposition rate, 
but we can be confident that this break is not caused by the 
increased count rate and insufficient pileup correction. 


3. DISCUSSION 
3.1. Soft X-Ray Emission 

Haberl et al. (2006) predicted a flattening of the 0.5-2. 0 keV 
light curve around day 7000 and predicted the light curve to 
reach its maximum around day 8000. Also, Borkowski et al. 
(1997) predicted the soft X-ray light curve to reach a maximum 
around day 9000, but with a flux of at up to a factor of 10 above 
the current X-ray flux. Indeed, Park et al. (2011) showed that 
the flux in the 0. 5-2.0 keV energy band flattened around day 
8200 and argued that this flattening is a result of the transmitted 
shocks entering a lower-density region. 

The improved calibration of the optical blocking filter now 
available removes a substantial part of this flattening. However, 
there still appears to be a break in the Chandra light curves, 
both in the piledup data (bare ACIS and zeroth-order HETG) 
and in the pileup-free dispersed data. Note that all the Chandra 
data presented here rely on the same calibration model for the 
optical blocking filter. The XMM-Newton EPIC pn data show no 
evidence for a sudden break around day 8200. Given that this 
break is visible in both dispersed and zeroth order, but not in the 
XMM-Newton EPIC pn data, we conclude that this is not likely 
an effect of an insufficient pileup correction, but rather is caused 
by a miscalibration of the increase in contamination on the 
optical blocking filters, combined with a change in instrument 
setup. Note that this contamination has a wavelength-dependent 
effect on the spectrum and hence has influence on the best- 
fit spectral parameters. Also, this contamination has less of an 
effect on the 3. 0-8.0 keV flux. 

3.2. Imaging 

Figure 1 shows the evolution of SN 1987 A’s morphology. To 
interpret these images, one has to keep in mind that some of the 
structure in the images might be an artifact of the deconvolution 
method we used to correct for the extension of the point-spread 
function of Chandra. Moreover, the point-spread function of 
Chandra is not entirely symmetric, as since 2000 December, 
there appeared a hook-shaped extension to the point-spread 
function, 111 which is not accounted for in MARX 4.4. This hook 
contains about 6% of the total flux and will probably show up 
after the deconvolution as a small extending structure to the 
remnant. 

The individual spots on the west side of the remnant appear to 
have merged in the last three epochs. One question that comes 
to mind when looking at Figure 2 is “how does this relate to 
the inner-ring structure as observed by Jakobsen et al. (1991)?” 
Ng et al. (2009) compared radii of X-ray and radio images of 
SN 1987 A and found comparable sizes, with radii larger than 
the radius of the equatorial ring in 2008. However, Dewey et al. 
(2012) showed a strong break in the evolution of the soft X-ray 
light curve after the shocks have overcome the equatorial ring, 
and we do not yet see such evolution in the 0. 5-2.0 keV light 
curve presented in Figure 4. 

To estimate how far the shocks progressed through the 
equatorial ring, we need to know when the interaction started. 
Sugerman et al. (2002) carried out a study on multiple epochs of 
the Hubble Space Telescope ( HST) images, ranging from 1994 
February and 2000 February. The first spot lit up by the blast 
wave is located in the northeast quadrant of the remnant and 
was detected in an HST observations taken 2933 days after 


16 See http://cxc.harvard.edU/ciao4.4/caveats/psf_artifact.html. 
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the explosion. Remarkably, this spot is located Cy/06 inside 
of the main ring as visible in optical (Sugerman et al. 2002). 
The second spot was only found in an HST observation taken 
1350 days later and by day 4999, 12 spots were visible in the HST 
observations (Sugerman et al. 2002). A lower limit to the width 
of the ring, as visible in the HST image, is of the order of O'.' 12 
(~10 12 km; Plait et al. 1995). Fitting a broken linear function 
to the radii of Figure 2 results in a velocity of 8500^ 2 j 3 9 oo km s' 1 
until day 5900^ 2 * * 2 8 3 q, after which the velocity is 1820 * 3 5 8 q km s _1 . 
Assuming that the shock wave reached the ring at day 4999 when 
12 spots were lit up (hereby ignoring the spots that were hit first, 
as those were protruding inward), the radius has increased by 
(1.17 ± 0.06) x 10 12 km (~120% of the width reported by Plait 
et al. 1995) since the interaction started. We used a Monte Carlo 
simulation to estimate the uncertainties on this value. Note that 
this distance relies on the assumption that the 12 spots reported 
by Sugerman et al. (2002) are part of the ring reported by Plait 
et al. (1995). 

Dwek et al. (2010) argued, based on the evolution of the 
infrared emission, that the forward shock only fully started to 
interact with the equatorial ring at day ~5600. This appears to 
be consistent with the increase in soft X-ray (lux around day 
6000 (Park et al. 2005). The latter was considered by Dewey 
et al. (2012) to be consistent with the lighting up of the opti- 
cal spots ~ 1 000 days earlier, as their models show some delay 
between the shock impact of dense material and the onset of 
X-ray emission. Nevertheless, if we take day 5600 as conser- 
vative assumption for the start of the explosion to interact with 
the ring, we estimate the increase of the radius of the X-ray 
emission since then to be 7.36 ± 0.05 x 10 11 km (~74% of the 
width reported by Plait et al. 1995). 

It is not straightforward to interpret the increase of the radius 
of the X-ray emission in terms of the three components as 
described in the introduction. The velocity is too high for the 
component with low (~0.3-0.5 keV) temperature (Zhekov et al. 
2009). Dewey et al. (2012) reported the contribution of a high 
velocity component to be ~20% in the HETG spectra, but 
its contribution to the broadband flux might be higher. Note 
that they measured the spectral lines to be much broader than 
1700 km s -1 . Therefore, the measured increase of the radius 
is likely a combination of the increase in radius of all three 
components. The radius of the X-ray emission is identified in 
Dewey et al. (2012) as following the progress of the forward 
shock/contact discontinuity as it moves through the uniform- 
with-clumps equatorial ring. Figure 9 (top) of Dewey et al. 
(2012) shows the forward shock/contact discontinuity locations 
along with the Racusin et al. (2009) radii and there is reasonable 
agreement. The speed of the forward shock there is of the order 
of 1700 knis" 1 . 

Current and future HST observations will reveal how far the 
shock has progressed, and how many of the optical spots have 
merged by now. A combined study with optical, X-ray, and 
infrared data will reveal a detailed picture of the structure of 
this equatorial ring, leading to better constraints on the last 
stages in the life of the progenitor of SN 1987 A. 

4. CONCLUSIONS 

We report our results for the three newest epochs of Chandra 

data on SN 1987A. Together with our reanalysis of older epochs, 

we paint a picture of the evolution of SN 1987A as witnessed by 

Chandra in the past 10 years. Based on this analysis, we reach 

the following conclusions. 


1 . The increase of the radius of the ring visible in X-rays since 
the explosion first interacted with the preexisting equatorial 
ring is 74%-120% of the lower limit to the thickness of 
this preexisting equatorial ring as measured by Plait et al. 
(1995). 

2. Given that the shocks have now traveled through a sub- 
stantial amount of the preexisting equatorial ring, a multi- 
wavelength study including high spatial resolution optical 
images is necessary to investigate how far the shock has 
traveled through the ring. 

3. The sudden break in the 0. 5-2.0 keV light curve around 
day 8200 reported by Park et al. (2011) was probably an 
instrumental effect. The absence of a break in the light 
curve likely also indicates that the shocks traveling through 
the dense parts of the ring have not yet fully overcome the 
extent of the ring. 
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